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Abstract 

In this paper we propose a new lifetime model, called the odd generalized exponential 
gompertz distribution, We obtained some of its mathematical properties. Some structural 
properties of the new distribution are studied. The method of maximum likelihood is used 
for estimating the model parameters and the observed Fisher’s information matrix is de¬ 
rived. We illustrate the usefulness of the proposed model by applications to real data. 
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1 Introduction 

In the analysis of lifetime data we can use the Gompertz, exponential and generalized exponential dis¬ 
tributions. It is known that the exponential distribution have only constant hazard rate function where 
as Gompertz, and generalized exponential distributions can have only monotone (increasing in case of 
Gompertz and increasing or decreasing in case of generalized exponential distribution) hazard rate. These 
distributions are used for modelling the lifetimes of components of physical systems and the organisms 
of biological populations. The Gompertz distribution has received considerable attention from demog¬ 
raphers and actuaries. Pollard and Valkovics fT4ll were the first to study the Gompertz distribution, they 
both defined fhe momenf generating funcfion of fhe Gomperfz disfribufion in ferms of fhe incomplefe 
or complefe gamma funcfion and fheir resulfs are eifher approximafe or leff in an infegral form. Lafer, 
Marshall and Olkin ifTOl described fhe negative Gomperfz disfribufion; a Gomperfz disfribufion wifh a 
negative rafe of aging paramefer. 

Recenfly, a generalization of fhe Gomperfz disfribufion based on fhe idea given in O was proposed 
by 0 this new distribution is known as generalized Gompertz (GG) distribution which includes the 
exponential (E), generalized exponential (GE), and Gompertz (G) distributions. A new generalization 
of th Gompertz (G) distribution which results of the application of the Gompertz distribution to the 
Beta generator proposed by [4], called the Beta-Gompertz (BG) distribution which introduced by ll^ . 
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On the other hand the two-parameter exponentiated exponential or generalized exponential distribution 
(GE) introduced by Q. This distribution is a particular member of the exponentiated Weibull (EW) 
distribution introduced by ifTll . The GE distribution is a right skewed unimodal distribution, the density 
function and hazard function of the exponentiated exponential distribution are quite similar to the density 
function and hazard function of the Gamma distribution. Its applications have been wide-spread as model 
to power system equipment, rainfall data, software reliability and analysis of animal behavior. 

Recently ifTSll proposed a new class of univariate distributions called the odd generalized exponential 
(OGE) family and studied each of the OGE-Weibull (OGE-W) distribution, the OGE-Erechet (OGE-Er) 
distribution and the OGE-Normal (OGE-N) distribution. This method is flexible because of the hazard 
rate shapes could be increasing, decreasing, bathtub and upside down bathtub. 

In this article we present a new distribution from the exponentiated exponential distribution and gompertz 
distribution called the Odd Generalized Exponential-Gompertz (OGE-G) distribution using new family 
of univariate distributions proposed by ifTSl . 

A random variable X is said to have generalized exponential (GE) distribution with parameters a, (3 if 
the cumulative distribution function (CDE) is given by 

F{x) = (1 - , X > 0, a > 0, /3 > 0. (1.1) 


The Odd Generalized Exponential family by ifTSl is defined as follows. Eet G(x; e) is the CDE of any 
distribution depends on parameter e and thus the survival function is G(x, e) = 1 — G(x; e), then the 


CDE of OGE-family is defined by replacing x in CDE of GE in equation (1.1) by 


Gix-,e) 

G{x,e) 


to get 


F{x-,a,e,l3) = 


_ £(^ 
1 _ g G(rr.e) 


X > 0, a > 0, e > 0, /3 > 0. 


(1.2) 


This paper is outlined as follows. In Section 2, we define the cumulative distribution function, density 
function, reliability function and hazard function of the Odd Generalized exponential-Gompertz (OGE- 
G) distribution. In Section 3, we introduce the statistical properties include, the quantile function, the 
mode, the median and the moments. Section 4 discusses the distribution of the order statistics for (OGE- 
G) distribution. Moreover, maximum likelihood estimation of the parameters is determined in Section 5. 
Einally, an application of OGE-G using a real data set is presented in Section 6. 


2 The OGE-G Distribution 


2.1 OGE-G specifications 

In this section we define new four parameters distribution called Odd Generalized Exponential-Gompertz 
distribution with parameters a, X,c,/3 written as OGE-G(0), where the vector 0 is defined by 0 = 

{a,X,c,l3). 

A random variable X is said to have OGE-G with parameters a,X,c,/3 if its cumulative distribution 
function given as follows 


F(x,e) 


i 


X > 0,a, X,c, (3 > 0, 


(2.1) 


where a, X, c are scale parameters and (3 is shape parameter. 
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2.2 PDF and hazard rate 

If a random variable X has CDF in (|2.1|) then the corresponding probability density function is 




f{x] 0) = a/3Ae“ec 


where x > 0, a, A, c, /3 > 0. 

A random variable X ~ OGE — G(0) has survival function in the form 


The failure rate function of OGE — G(0) is given by 




h{x) = 


fix) 

Six) 




^-1 


-ale 

1 — e L 




0-1 


1- ll-e~ 




( 2 . 2 ) 


(2.3) 


(2.4) 




Figure 1: The CDF of various OGE-G distributions for some values of the parameters. 




Figure 2: The pdf’s of various OGE-G distributions for some values of the parameters. 
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Figure 3: The hazard of various OGE-G distributions for some values of the parameters. 


3 Statistical properties 

In this section, we study some statistical properties of OGE-G, especially quantile, median, mode and 
moments. 


3.1 Quantile and median of OGE-G 

The quantile of OGE-G(0) distribution is given by using 

F{xq) = q. 


Substituting from (2.1) into (3.11, can be obtained as 


Xo = - In 1 + — In 


1 


1-ln(l — qP) 


a 


, 0 < g < 1. 


Setting g = ^ in (3.2 1 , we obtain the median of OGE-G(q!, A, c, /3) distribution as follows 


1 f c 

Med = - In < 1 + — In 
c A 


a Z 


(3.1) 


(3.2) 


(3.3) 


3.2 The mode of OGE-G 

In this subsection, we will derive the mode of the OGE-G(0) distribution by deriving its pdf with respect 
to X and equal it to zero thus the mode of the OGE-G(0) distribution can be obtained as a nonnegative 
solution of the following nonlinear equation 


c + Ae“-aAe“e^("“-i) 


-a e 

1 — e L 

A f^cx _ 


'-1 


+ 




aA(/3 — l)e‘^*e 


0 . 


(3.4) 


It is not possible to get an explicit solution of (3.4i in the general case. Numerical methods should be 
used such as bisection or fixed-point method to solve it. 
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3.3 The moments 


Moments are necessary and important in any statistical analysis, especially in applications. It can be 
used to study the most important features and characteristics of a distribution (e.g., tendency, dispersion, 
skewness and kurtosis). In this subsection, we will derive the rth moments of the OGE-G(0) distribution 
as infinite series expansion. 

Theorem 3.1. If X OGE — G(0), where 0 = (a, A, c, /5), then the rth moment of X is given by 
' = ^ (_-\ f+j+k+m f + lyjj - k + l)^r! 

G 1 • J \rn) [kj c^+r+ijiiU£_rn + lY+^ ' 

i=0 j=0 k=0 e=0 m=0 ^ \ \ j \ ' J 

Proof The rth moment of the random variable X with pdf f{x) is defined by 


= / x'^f{x)dx. 

Jo 


Substituting from (2.21 into (3.51, we get 






/3-1 


dx. 


(3.5) 


(3.6) 


Since 0 < 1 — e 


—a lee'' 


< 1 for X > 0, we have 

/3—1 oo 




(3.7) 


Substituting from (|3.7[) into p.6|), we obtain 


tf', = p^apx{-iy 


Using the series expansion of e ^ we get 

' a\f 1 \i+j f~ of {i + 1 )-^ 

^ = EE“W-iU( i 

i=0 j=o ^ ^ 

Using binomial expansion of (e c — 1)7, we get 

oo oo i 

br = 

i=0 1=0 k=0 

Using series expansion of ec vve get 

OO OO J OO 

Mr = 




X e ec 


— ly dx. 


4=('’ 7') (7) ?fEm±dl r 

«_n -v—n G—n \ / \ / J' */0 


prgCXgAQ_fc+l)(e.._l)^X. 


EEEE(-iU’*‘ ('’■') (i) " 

i=0 j=Q k={) 1=^ ^ / \ / J' ' 

roo 

/ x"e'^(e'^ - l)^dx, 

E E E E E (- y 7 0 (™) ( 7 ) 

„_n A—n u—c\ D—Ci -v>o—n \ / \ / \ / J* • 


2=0 j=0 /c=0 £=0 m=0 
roo 

Jo 
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By using the definition of gamma funetion in the form 


r( 2 ;) = ^dt, z,x > 0 


Thus we obtain the moment of OGE-G(0 ) as follows 


oo cx j oo £ / a -I 

^_-^y+j+k+m fP 1 

i=0 j=0 k=0 £=0 m=0 ' 

This completes the proof. 


£\ f j\ + l)^{j — k + l)^r! 


mj \kj c^+^+^j!.^!(£ — m + 1)^+^ 


□ 


4 Order Statistic 

Let Xi,X 2 , • • • , Xn be a simple random sample of size n from OGE-G(0) distribution with cumulative 
distribution function F(x; 0) and probability density function /(x; 0) given by (2.1 1 and (2.2 1 respec¬ 
tively. Let < X 2 :n < ''' < Xnm denote the order statistics obtained from this sample. The 
probability density function of Xr-.n is given by 


fr-.n{Xj 0) — 


1 


[F(x,0)rMl-F(x,0r-'’/(x;0), 


B{r,n — r + 1) 


(4.1) 


where /(x; 0) and F(x,0) are the pdf and cdf of OGE-G(0) distribution given by (2.1 1 and (2.2 1 
respectively and B {.,.) is the beta function. Since 0 < F{x, 0) < 1 for x > 0, we can use the binomial 
expansion of [1 — F{x, 0)]”~^ given as follows 


[1 - F(i, e)]"-' = y; 


2=0 


n — r 


i-iy [F{x,e)Y 


(4.2) 


Substituting from (4.21 into ( |4.1| ), we have 

fr:n{Xj 0) — 


B{r, n — r + 1) 


/(x;0)^ 


2=0 


n — r 


i-iy [F(x,0)] 


2 +r—1 


(4.3) 


Substituting from (2.11 and (2.2 1 into (4.31, we obtain 

(—l)*n! 


fr-.nix; a, A, c,P) = Y, 
i=0 


— l)!(n — r — iy.{r + i) 


-f{x,a,\,c, (r+ i)/3) 


(4.4) 


Thus fr-.n{x] a, A, c, /3) defined in (4.4) is the weighted average of the OGE-G distribution with different 
shape parameters. 


5 Estimation and Inference 

Now, we determine the maximum-likelihood estimators (MLE’s) of the OGE-G parameters. 
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5.1 Maximum likelihood estimators 

Let Xi,X 2 , - ■ ■ ,Xn be a random sample of size n from OGE-G(0), where 0 = (a, A,c, /3), then the 
likelihood function £ of this sample is dehned as 




(5.1) 


2=1 


Substituting from (2.21 into (4.31, we get 


C = Y[ a(3\e^ 


D _i] | ^ [e c D-l]^ 


2=1 


The log-likelihood function is given as follow 


L = nln(a) + nln(/3) + nln(A) + c^Xj H—— 1) — a 


e c 




-1) 


1=0 


i=0 


i=0 




i=0 


(5.2) 


The log-likelihood can be maximized either directly or by solving the nonlinear likelihood equations 


obtained by differentiating equation (5.2i with respect to a, A, c and (3. The components of the score 
vector [/(0) = (f^, If, f,§) are given by 


OL 

da 

dc 


n 




n 

a 


i=\ 

n 




(5.3) 


2=1 


A(e-i-l) 


l)-(/3-l)^^ 


1 


2=1 


^-1 


(5.4) 


- 1 


n 1 


A c 


2=1 


+ i _ 1) _ “ ~ 


ai/3 - 1) A - 1) 


2=1 


1=1 e L J — 


(5.5) 


^ Xj - ^ T{xi, A, c) - q; ^ r(xi, A, c)e = 




+ (/5 - 


r(xi, A, c)e = 


2=1 2=1 


2 = 1 


2=1 e L 


- 1 
(5.6) 


where r(xi, A, c) = - 1) + ^Xje^L 

The normal equations can be obtained by setting the above non-linear equations (5.31-( 5.6 1 to zero. That 
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is, the normal equations take the following form 


n 


!>= 0 , 


2=1 


n n . 

n , A 


a 




He‘^^i-1) _-l] _ (R _ 


ec 

2 = 1 

1 ^ 

n 1 a 




1 _ 




= 0 , 


(5.7) 

(5.8) 


A c 


+ - 5^(e“‘ - 1) - - J^(e“‘ - + 


2=1 


2 = 1 


a(/3- 1) A (e^^* - 

*=i e L J — 1 


= 0 , 


(5.9) 

n n n ^ i \ \ -(e‘^^i-1) 

^ ^ r(xi, A, c) - a ^ r(xj, A, c)ec + (/S - !)« ^ ^ -;-= 0. 


2 = 1 2=1 


2=1 


1=1 e L J — 1 


(5.10) 

The normal equations do not have explicit solutions and they have to be obtained numerically. From 


equation (5.71 the MLEs of /3 can be obtained as follows 

(!= - 


Et., li 1 - e 


—a I e c ^ 


(5.11) 


Substituting from ( |5.11| ) into (5.8 1 , (5.91, and (5.101, we get the MLEs of a, A, c by solving the following 
system of non-linear equations 


n . n 


a 


i=l 


i=l « 
e L 




(5.12) 


n 1 

“ + T 

A c 


2=1 


+ 1 - 1 ) - “ + 


J - 1 

(5(7 — 1) ^ (e'^^' — l)e^ 


E 




i=l 


i=l “ 
e L 


'-1 


= 0 , 


J - 1 


(5.13) 


^ ^ T{xi, A, c) - (5 ^ r(xi, A, 


c ec 




0 - l)aX] 


T(xi, A,c)e7'=‘'''‘ 1) 


2=1 2=1 


2=1 


2=1 » 
e 


Xf^cx. 

e c 




= 0 , 


J - 1 


(5.14) 


where r(xj, A, c) = “ 1) + . 

These equations cannot be solved analytically and statistical software can be used to solve the equations 
numerically. We can use iterative techniques such as Newton Raphson type algorithm to obtain the 
estimate $. 

5.2 Asymptotic confidence bounds 

In this subsection, we derive the asymptotic confidence intervals of the unknown parameters a, A, c, /3 
when a>0,A>0, c>0 and /? > 0. The simplest large sample approach is to assume that the 
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MLEs(q!, X,c, /3) are approximately multivariate normal with mean (a, X,c, /3) and eovarianee matrix 
where is the inverse of the observed information matrix whieh defined as follows 


4"' = - 


92 ^ 

d^L 

d^L 

d^L 1 

- 

1 

~dcA 

dXda 

dcda 

dpda 



d^L 

d^L 

d'^L 

d^L 



dadX 


dcdX 

dfSdX 



d^L 

d^L 

d^L 

92 



dadc 

dXdc 

IkA 

dpdc 



92 

92 

d^L 

92 



dad/S 

dXdl3 

dcdl3 

9/32 



var{Q:) 

COV 

(A, a) 

cov{c, d) 

cov{j3, d) 

cov{a, X) 

var(X) 

cov{c, / 


cov{(3, A) 

cov{d,c) 

COV 

(A,c) 

var{c 


cov{l3, c) 

cov{d, f3) 

COV 

(4/3) 

cov{c, / 

3) 

var{j3) 


The seeond partial derivatives ineluded in /g ^ are given as follows 
d^L _ —n 

d^L _ A )-lje L J 

dadl3 ^ 

1=1 1 _ g L J 

d^L _ a A - 1) 

dXdp c ^ ’ 

1=1 1 — e L J 

52 L A AtBt 


r A(eCa!i 
—a e c '• 


—a e c — 1 


f, 1 - ’ 

ta J 




cla2 q,2 


(<9-l)E 


—a ^>—1 

— el -I 


52 L 

dXda 


dcda 


= — ^ - 1) e = + — _ — y~] 






A-iBiCi 


Afe^^ —li 1 ' 
—a ^''-1 

— el J 


(5.15) 


U ij —Tl CX ,9 

2=1 


2 A(eCa:i_l) Q!(/3 — 1) — 1) AiDi 


:r|i_^-„[e4--)-r]y 


d'^L 1 


= -xT.b,--X + + 

2=1 2=1 ^ 

" s jlI 


- 1 ) 


I (/3 - l)a y-^ 


2 
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d^L 

dc^ 




i=l 


i=l 




Bj Di Ei I 1 — 




A 


2 = 1 


Ji V 


where ^-1, i?* = — 4(6“^* — 1) + 


Ct = l- 




Dt = 1 — ae 






Ei = ^ - 1 ) - 

C'^ c 


The asymptotic (1—7)100% confidence intervals of a, A, c and (3 are a^nz-i y/var{a), A±Z 2 ^ var{X), c± 
2;2 y/var{c) and ^ ± zi \Jvar0)) respectively, where zz is the upper ( ^)th percentile of the standard 


2 ' 2 

normal distribution. 


6 Data Analysis 

In this section we perform an application to real data to illustrate that the OGE-G can be a good lifetime 
model, comparing with many known distributions such as the Exponential, Generalized Exponential, 
Gompertz, Generalized Gompertz and Beta Gompertz distributions (ED, GE, G, GG, BG), see ||3l|6l|7llll. 
Consider the data have been obtained from Aarset |T|, and widely reported in some literatures, see 
ll2ll3l|9l[l2j[T3l. It represents the lifetimes of 50 devices, and also, possess a bathtub-shaped failure rate 
property. Table 1. 


Table 1: The data from Aarset ifTl. 


0.1 

0.2 

1 

1 

1 

1 

1 

2 

3 

6 

7 

11 

12 

18 

18 

18 

18 

18 

21 

32 

36 

40 

45 

46 

47 

50 

55 

60 

63 

63 

67 

67 

67 

67 

72 

75 

79 

82 

82 

83 

84 

84 

84 

85 

85 

85 

85 

85 

86 

86 



Based on some goodness-of-fit measures, the performance of the OGE-G distribution is compared with 
others five distributions: E, GE, G, GG, and BG distributions. The MEE’s of the unknown parameters for 
these distributions are given in Tables 2 and 3. Also, the values of the log-likelihood functions (-L), the 
statistics K-S (Kolmogorov-Smirnov), AIC (Akaike Information Criterion), the statistics AICC (Akaike 
Information Citerion with correction) and BIG (Bayesian Information Criterion) are calculated for the 
six distributions in order to verify which distribution fits better to these data. 
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Table 2: MLE’s of parameters for Aarset data. 


Model 

a 

$ 

A 

c 

-E 

E 

0.0219 

- 

- 

- 

241.0896 

GE 

0.0212 

0.9012 

- 

- 

240.3855 

G 

- 

- 

0.00970 

0.0203 

235.3308 

GG 

- 

0.2625 

0.00010 

0.0828 

222.2441 

BG 

0.2158 

0.2467 

0.00030 

0.0882 

220.6714 

OGE-G 

0.0400 

0.1940 

0.000345 

0.0780 

215.9735 


Table 3: The AIC, CAIC, BIC and K-S values for Aarset data. 


Model 

AIC 

AICC 

BIC 

K-S 

p-value(K-S) 

E 

484.1792 

484.2625 

486.0912 

0.19110 

0.0519 

GE 

484.7710 

485.0264 

488.5951 

0.19400 

0.0514 

G 

474.6617 

475.1834 

482.3977 

0.16960 

0.1123 

GG 

450.4881 

451.0099 

456.2242 

0.14090 

0.2739 

BG 

449.3437 

450.2326 

456.9918 

0.13220 

0.3456 

OGE-G 

423.9470 

424.8359 

447.5951 

0.13205 

0.3476 


Based on Tables 2 and 3, it is shown that OGE-G(a, A, c, /3) model is the best among of those distribu¬ 
tions beeause it has the smallest value of (K-S), AIC, CAIC and BIC test. 

Substituting the MEE’s of the unknown parameters a, A, c, /3 into (5.151, we get estimation of the vari¬ 
ance covariance matrix as the following 



1.5022 X 10-3 

2.7638 X 10-7 

-1.4405 X 10-^ 

9.6331 X 10-^ 

Io^ = 

2.7638 X 10-^ 

4.4773 X 10-® 

-1.7837 X 10-6 

2.785 8 X 10-6 

-1.4405 X lO-'^ 

-1.7837 X 10-6 

8.4428 X 10-6 

-1.9127 X 10-^ 


9.6331 X lO-'^ 

2.7858 X 10-6 

-1.9127 X 10-^ 

1.5309 X 10-3 


The approximate 95% two sided confidence intervals of the unknown parameters a, X,c and /3 are 
[0,0.116], [0,0.001], [0.060,0.096], [0.116,0.270], respectively. 

To show that the likelihood equation have unique solution, we plot the profiles of fhe log-likelihood 
function of a, A, (3 and c in Eigures 4-5. 




Eigure 4: The profile of fhe log-likelihood funclion of a, A. 
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Figure 5: The profile of the log-likelihood function of c, (3. 

The nonparametric estimate of the survival function using the Kaplan-Meier method and its fitted para¬ 
metric estimations when the distribution is assumed to be ED, GED, GD, GGD and OGE — GD are 
computed and plotted in Figure 6. 



Figure 6: The Kaplan-Meier estimate of the survival function. 

Figure 7, gives the form of the hazard rate for the ED, GED, GD, GGD, BGD and OGE — GD which 
are used to fit the data after replacing the unknown parameters included in each distribution by their MLE. 



Figure 7: The failure rate function for the data. 
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7 Conclusions 

In this paper, we propose a new model, ealled the Odd Generalized Exponential-Gompertz (OGE-G) 
distribution and studied its different properties. Some statistieal properties of this distribution have been 
derived and diseussed. The quantile, median, mode and moments of OGE-G are derived in elosed forms. 
The distribution of the order statisties are diseussed. Both point and asymptotie eonfidenee interval 
estimates of the parameters are derived using the maximum likelihood method and we obtained the 
observed Eisher information matrix. We use applieation on set of real data to eompare the OGE-G 
with other known distributions sueh as Exponential (E), Generalized Exponential (GE), Gompertz (G), 
Generalized Gompertz (GG) and Beta-Gompertz (BG). Applieations on set of real data showed that the 
OGE-G is the best distribution for fitting these data sets eompared with other distributions eonsidered in 
this artiele. 
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